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ABSTRACT 


Part A of this report examines several high-frequency models for 
nonprincipal-plane scattering from a rectangular, perfectly conducting 
plate. Two methods, the Method of Equivalent Currents and corner 
diffraction coefficients, are considered. Formulations for 
second-order Physical Theory of Diffraction equivalent currents and 
for corner diffracted fields are presented. Comparisons are made 
among the following plate models: first-order Physical Optics 
equivalent currents, first-order Geometrical Theory of Diffraction 
equivalent currents, first-order Physical Optics/Physical Theory of 
Diffraction equivalent currents, second-order Physical Theory of 
Diffraction equivalent currents, corner diffraction coefficients, 
Moment Method, and experimental results. Results away from grazing 
are accurate using only first-order terms. Near grazing, second-order 
and corner diffraction terms improve the results for many cases . 

Part B of the report investigates the pattern control of horn 
antennas using lossy materials to coat the inner walls of the horn. 
Integral Equation and Moment Method techniques are used to formulate 
the problem. It is clearly demonstrated that side lobe level 
reduction can be achieved using impedance surfaces on the inner walls 
of the horn. 
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In this report, two new models of the plate for 
nonprincipal-plane scattering are explored. The first is a revised 
version of the PO/PTD model with second-order PTD equivalent currents 
[6] included to account for second-order interactions among the plate 
edges. The second model uses a heuristically derived corner 
diffraction coefficient [7], [8] to account for the corner scattering 
mechanism. The patterns obtained using the newer models are 
comparedto the data of the previously reported models, MM, and 
experimental results . 
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II. Theory 


A. Second-Order PTD Equivalent Currents 

*s „ith many versions of the MEC, the second-order PTD equivalent 

currents are formulated using the canonical perfectly conducting, 

, • w used to approximate other geometries, 

infinite wedge geometry which is used PP 

For the rectangular plate shown in Fig. 1, each edge is modeled as 
infinite half-plane by setting the exterior wedge angle to 2n. 
i, a valid approximation as long as the edges are electrically 
isolated; thus, the accuracy of this model increases as the electrical 
Sire of the plate increases. The general wedge geometry is shown in 
Fig. 2. This geometry is applicable to both first- and second-order 
equivalent currents. The directional vectors and angles are; 

S' = the unit vector in the direction of incidence. 

I = the unit vector in the direction of observation. 

t = the unit vector tangent to the edge of 

directed so that it encircles the sc.tterer in 

counterclockwise manner. 

yv pdae of interest, lying 

n = the unit vector normal to the eage 

on the illuminated face. 

A 

fS ' = the angle between s' and the edge. 

A 

(3 = the angle between s and the edge. 

2 = the angle between the illuminated face and the 

edge-fixed plane of incidence. 

# . the angle between the illuminated face and the 

edge-fixed plane of observation. 
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a — the skew angle of integration across the surface. 


In terms of the directional vectors, the directional angles 

-l A A 

/3 0 ' = cos (s' • t) 

_ ^ A A 

|3 0 = COS (S • t) 


0 ' 



COS 


s'X t 

A A 

| s'X 1 1 



: 


are : 

(la) 

(lb) 

(lc) 


-1 s X t 
0 - COS — ~ 

U |» X t| 

For monostatic scattering, ip * 
flat plate, the unit vector 
geometry considered is: 



ip' f 13 Q - tt - 0 O ' , and s - -s' . 
in the direction of incidence 


(Id) 

For the 
for the 


A A 

s ' = -a r 

A A 

= -a x sin0' cos 0' - a y sin0'sin0' - 

A 

a z cos0' 

(2) 

A different 

set of directional vectors and 

angles 

must be 

formulated for each edge. The vectors for each edge 

are : 


Edge 1 : 

A A A A 

t x = a y n x = -a x 


(3a) 

Edge 2 : 

A A A A 

t 2 = -a y n 2 = a x 


(3b) 

Edge 3 : 

A A A A 

t 3 = -a x n 3 = -a y 


(3c) 

Edge 4 : 

A A A A 

t A = a x n 4 “ a y 


(3d) 

The resulting (3 Q f 

functions for each edge are: 



Edge 1 : 

cos/8 0 ^ = -sin0'sin$ r 


(4a) 


sin0 o ' = J 1 - sin 2 0'sin 2 0' 


(4b) 

Edge 2 : 

COSp D ^ = -COSp D ' 


(4c) 


sin/3 0 ^, = sinP D ^ 


( 4d) 
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Edge 3: cos/3 Q ' = sin0'cos<£' 


/ 2 2 
1 - sin 0' cos <p f 


sin ^ 


Edge 4: cos£ Q ' = -cosf3 0 ' 

sin £ = sin P 0 3 


Using the /3 0 ' definitions, the 0' functions can be expressed as: 

. . -sin0 f cos0' 

Edge 1: cosi = 


sin 0 1 = 


sin£ 0 ^ 

COS0' 


sin/3 0 ' i 

Edge 2 : co s0 2 “ “cos^ 

sini// 2 " sini// 1 

_ -sin0'sin$' 


Edge 3 : 


00501 ■ =in|3 0 ; 


sin ^ 3 = 


COS0' 

sin 0 o . 


Edge 4 : 


(4e) 
( 4 f ) 
(4g) 
(4h) 

(5a) 

(5b) 

(5c) 
( 5d) 

(5e) 
( 5 f ) 


( 5g) 
(5h) 


cosi/i 4 = -cos^ 3 
sin^ 4 = sini// 3 

For each edge the tangential components of the incident electric 

and magnetic fields are needed to determine the corresponding 

equivalent currents. Both soft and hard polarizations are considered. 

The incident fields are: 

Soft Polarization 
_ i ~ -jk*r 

- = a (f> E ° ® “ (6a) 

H 1 = a Q (1 /tj) E 0 e' 3 -*- (6b) 
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(6c) 


Hard Pol^£il£tion 
- jk. • r 

E = a 0 E o e 


(6d) 


- jk*£ 

H 1 - (1/T,) E ° 6 


« “ . the incident fields 

To simplify the rectangular plate analyst 

d to the rectangular coordinate system. 
are transformed 

vector, £, is: ^ (7) 

r - a x x + a y y + a z z 

, nr k for the incident field is: 

The propagation vec a . + a cos0' ) (8) 

k( ; sin0' cos<p' + a y sin 0 'sin<^ a z 

* = -k a r = ■*< x incident fields 

ljr coordinate system, th 
to the rectangular cooru 
With respect to 


are : 


Soft Polariz ation 

in0’ coslp' +ysin0' sin</>' +zcos0 ) 


i jklxs 

E = E 0 e 


X l- a > 


A 

sin <p' +a y cos <f> 1 


i 1 = (1/7)) E 0 e 


jU xsin0'cos0'tysin0'sin<#.' + zcos 


9 ') 


x [a K cos0' cos*' +a y cos9' sin^' +a z sin0' 1 


(9a) 


(9b) 


Hard Polarization 

<sin 0'ccs0'tysin0'sin^ + zcos0 ) 


Ei - E 0 e 


jk (xs 


X 


[a cos 0 'cos^»'+a y 


C os0'sin4>'+a z sin0'l 


H i = - ( l/*n) E 0 e 


x [-a x sin^+a y cos0'] 


(9c) 


( 9d) 
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The scattered fields for far-field observations are expressed in 
terms of the vector potentials: 


E r 

= 0 




(10a) 

E e 

£ -jo> 

IA e + 1 



(10b) 

E 4> * - ^ F e J 

The vector potentials are: 



(10c) 

A = 

; J£ 

[ I(x',y',z') 

e- jkr 

d V 

(11a) 


47T 

J 

( 

R 

F = 

_£ | 

f 

M(x' ,y' , z' ) 

-jkr 

e 

di' 

(lib) 


4tt J 

R 


c 


For far field analysis, the following simplifications can be made: 

A 

R = r ~ r'cost? = r - r' • a r (for phase variation) (12a) 

R - r (for amplitude variation) (12b) 

where 

A /V 

r' = a x x + a y y ( 1 3) 

for the flat plate oriented as in Fig. 1. For phase variation: 

R = r - xsinG'cos#' - ysinG'sin^' ( 14 ) 

Finally, the vector potentials for far-field scattering from the 
plate become : 


A = 


M e 


- jkr 


4tt r 


T/ v jk (xsinG' cosd)' +ysinG' s in$' ) 

I <X, y, z) e r Y 


dt 


-jkr 


F - 


e e 
471 r 


| M(*,y, 


jk (xsinG' cos(b f +ysinG' sin<b' ) n 
z) e ^ 1 Y dt 


(15a) 


(15b) 
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The integrals are evaluated along the perimeter of the plate in a 


counterclockwise manner. On 

the surface of 

the plate. 

z=0 . Also, 

for each edge, the 

coordinates 

are : 



Edge 1 : 

x = a 

XI 

VI 

>i 

VI 

XI 

1 


(16a) 

Edge 2: 

x = -a 

-b £ y 3 b 


(16b) 

Edge 3 : 

-a < x ^ a 

y - b 


(16c) 

Edge 4 : 

-a < x < a 

y = -b 


( 1 6d) 

To simplify the 

derivation. 

each vector 

potential 

integral is 


represented as a sum of four integrals, each corresponding to an edge 
of the plate. For the second-order currents, numerical integration is 
used in evaluating the integrals . ][ and M must be determined 
separately for each edge taking into account the individual 
geometries . Again, each edge is viewed as the truncation of an 
infinite half-plane . 

The general geometry for the formulation of the second-order 
currents is shown in Fig. 3. The directional vectors and angles are 
for the edge of first-order diffraction defined above in the 
description of Fig. 2. The diffracted ray travels along cr at a skew 
angle of fS Q r . The axis of integration across the structure is again 

A 

<r, which effectively eliminates all singularities except the Ufimtsev 
singularity for forward observation at grazing incidence. The 
distance from the first point of diffraction, 0 lf to the second point, 
C> 2 , is if. The tangent vectors to the edge of first diffraction is t ^ 
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4 y 


/ 


4 y 



Fig. 3. Geometry for second-order components. 
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, ■ riiupop a j -0 oriented so that 

and the vector at the second edge is 2 - 

they encircle the scatterer in a counterclockwise manner. 

The second-order PTD equivalent currents components are 

defined as : 


f 

Mo = 


l 2 “ *2 

sin /3 2 


- 


2 

sin fi 2 


where 


[ s x K 2 ] 


[ (t, X S) • ] 


£ 


( 17 a) 


( 17 b) 


K 2 « |cr x t- 


-jk£ («••■> ( -f (<r * s) 


j e 


dcr 


( 17 c) 


The surface fringe current density, jt is expressed in terms of the 
exact fringe-current scattering solution to the wedge; therefore, K 2 
consists of a contour integral in the complex plane that is integrated 
along ff and evaluated only at the upper endpoint. K 2 is evaluated 
using the Method of Steepest Descent and the result for a half-plane 

is [5]: 

- jkje 


*2 = 


AV ~2 Iff X t; 


jksin 2 p o ' (n+cosi//' ) 


F [V L( 1 ~H) ] 
V 1-u 


X 


[nj (t x * H„)cos(0'/2) + t^t!- Hq) (icotp 0 ' cos (</*' /2) 

+ tjlti-Ej) (1/7)) ( 1 -M>csc 0 o 'sin( 0 ' /2) ] 


_ F [V L(l + COS 0 r T Hq) - t x <t x * H„) COt/3 0 ' COS0' 
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(18) 


+ titt-* Eq) ( 1 /tj) csc|3 0 ' sin0' ] 


The expressions 


for the currents simplify to: 


f * f sin \p 2 

M 2 = -D n 2 * K2 


(19a) 


I* - t 2 > kJ - n 2 * £2 cotp 2 cosi/f 2 


(19b) 


where K 2 is 


f • given in (18) for a half-plane. Other necessary 


quantities are: 


A 



1 2 


n l' 


n 2 


the unit vectors tangent to the edges of first 
and second-order diffractions, respective y, 
oriented so that they encircle the scatterer in a 
counterclockwise manner. 

the unit vectors normal to the edges of first and 
second-order diffractions, respectively, pointing 
inward and lying on the illuminated surface of the 
scatterer . 


A 

cr 


the unit vector in the direction of integration 
at an angle 0 O ' with respect to the edge of the 
so that it represents the grazing diffracted ray. 


skewed 

wedge 


or 


or 


* ~ 0 , (20a) 

cr = n x sin/3 0 ' + t x cos/3 0 ' 

jg = the distance along cr from O x to 0 2 , the points of 
first- and second-order diffraction, respectively. 


L = k£ sin 2 p o f 


F(x) = the modified Fresnel transition function. 


F (x) 


V j / tt e^ x 



x 


(20b) 


(20c) 
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cosy - COS P n ' _ 

M 2 

sin |3 0 ' 


sin ( 2 r/ 2 ) 


(20d) 


sin p 0 ' 

A A 

■y = the angle between cr and s. 

cosy = cr • s = 3infi 0 ' sinPjCOSi/^ + cos/3 0 ' cos^ 


(20e) 


0if 0, = the polar angles of s in the coordinate system 
local to the edge of interest, the n^y,!* 
coordinate system. 


where 


cos/^ = s • t i 

A A 

sin/^cos^i - s • n i 

A A 

sinfi i sin\p i = s • y 


( 2 0 f ) 
(20g) 
(20h) 


The directional angles 0 O ' , 3 0 , <P' , and 0 are defined in (4) and (5). 

The second-order diffracted field is obtained by substituting M 2 
and I f 2 from ( 19 ) into the vector potential integrals of (ID, 
integrating, and then substituting into (10) . Numerical integration 
must usually be used to find the integrals of (11) . The limits of 
integration on the integrals of (11) are found using ray tracing. The 
area that the first-order fields affect is bounded by the two extreme 
first-order diffracted rays. Fig. 4 illustrates this procedure. Edge 
AB illuminates the curve from A' to B' . Integration is along the 
boundary from A' to B' . Often illuminated regions overlap due to 
interactions from many edges. 

The total first- and second-order fields are found by adding the 
fields due to scattering by the PO currents, the PTD components, and 
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the second-order components. These fields are valid for all 
directions of illumination and observation, on and off the Keller cone 
of diffracted rays, except for the forward direction due to grazing 
incidence where an infinite singularity exists. This is the Ufimtsev 
singularity. 

In terms of the directional vectors, the directional angles are: 


are : 


A 

cos/3 2 = s 

* t 




(21a) 

sin/3 2 cosi// 2 

A 

= S * 

A 

n 2 



(21b) 

sinf^ 3 ^ 1 /^ 

A 

= s • 

A 

y 



(21c) 

edge of second- 

-order 

diffraction, the directional 

vectors 

Edge 1 : 

A 

t 21 = 

A 

a y 

A A 

n 21 ~ “ a x 


(22a) 

Edge 2 : 

A 

t 22 = 

A 

-a y 

A A 

n 22 = a x 


(22b) 

Edge 3 : 

A 

t 23 “ 

A 

_a x 

A A 

n 23 = _a y 


(22c) 

Edge 4 : 

A 

t 2 4 = 

A 

a x 

A /V 

n 24 = a y 


(22d) 

definitions of 

(22) 

in the equations of (21) 

along 

with the 




sin0 2 

cot/3 2 


n of s' 

from 

(2) , 

the factors sin/3 2 and 

COSl/^2 

in (19) 


simplify to the following for each edge: 

sin^i _ cos0' 


Edge 1 : 


(23a) 


sinp 21 sinVcosV+cosV 


-sin 0 r cos0'sin0' 

cot|3 21 cost p 21 - j 2 

sin 0'cos <f> r +cos 0' 


(23b) 
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Edge 2 : 


Edge 3 : 


Edge 4 : 


sin0 2 2 

cos0 f 

(23c) 

s in(3 22 sin 2 ! 

2 2 

9' COS (p f +COS 0' 


COt/3 2 2 COSI ^22 “ 

-sin 2 0' cos <p' sin <f>' 

(23d) 

2 2 2 

sin 0'cos 0'+cos 0' 

sini // 2 3 

COS0' 

(23e) 

sinp 23 . 2 

1 2i sin 

2 2 
0' sin <f>' +cos 0' 


c°t/3 2 3 cosl / , 2 3 = 

2 

sin 0' cos0' sin0' 

(23f ) 

2 2 2 
sin 0'sin <p r +cos 0 f 

sini/f 24 

COS0 f 

(23g) 


sin/3 24 


sin 2 0' sir\ 2 <p r +cos 0' 


2 

sin G' cos <t>' sin<p r 

COtP 24 COSl// 2 4 = 2 2 277 

sin 0' sin +cos 0 


(23h) 


The remaining factor in the equivalent currents equations that 
must be determined for each edge is K f 2 . The preceding vectors and 
functions refer to the edge of second-order diffraction only. The 
vectors and functions involved in the definition of involve the 


edge of first-order diffraction only except for the \<r X t 2 1 factor 
which involves directional vectors from the edges of both first- and 
second-order diffraction. For a half plane, k| is given in (18). 

The t x and n x vectors are defined for each edge in (3) . The 

incident fields, H 1 and E 1 are defined in (9). F(x) is the modified 
Fresnel transtion function of (20c) . The directional angles 0 O ' and 0 
are the same as those introduced at the beginning of this section. 
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All the necessary functions of these angles derive easily from the 


relations of (4) and 

(5) . 


The unit vector 

in the direction of integration, <r, for each 

edge 

is : 



Edge 1 : 

cr 1 = -a x sin 0 o ' + a y cosp 0 ' 

(24a) 

Edge 2 : 

o~ 2 = a x sin3 02 - a y cos/3 0 ' 

(24b) 

Edge 3 : 

a * A , 

cr 3 = -a x cos 0 o ' - a y sin/3 03 

(24c) 

Edge 4 : 

= a x cosp D ' + a y sin 0 o ' 

(24d) 

In general the 

term /i is defined as: 


cos y ‘ 

2 

- COS /3 0 ' 

(25) 

* 2 

sin 0 O ' 


where 



A 

COST “ O' 

a. 

• s 

(26) 

For each edge, the 

y functions are: 


Edge 1 : 

cos*! = -sin 0 o 'sin 0 'cos 0 '+cosp o 'sin 0 ;sin 0 ' 

(27a) 

Edge 2 : 

cosr 2 = sin 0 o ' sin 0 ' costf>' -cos 0 o ' sin 0 ' sin<p' 

(27b) 

Edge 3 : 

cosr 3 = -sin 0 o ' sin 0 ' cosf -cos 0 o ' 3 sin 0 ' sinf 

(27c) 

Edge 4 : 

cosy,) = Sin0 o 'sin0'cos^'+cos0 o 'sin0'sin0' 

(27d) 

Using these definitions and those of (4), the M's for each edge 

are : 



Edge 1 : 

2 

cos y ^ - cos /3 0 ^ 

Mi 5 

sin 0 O ' 

(28a) 

Edge 2 : 

2 

COS^ 2 " cos $o 2 

^2 2 

sin f3 0 ' 2 

(28b) 


-2 0 - 



Edge 3 : ^3 " 


Edge 4 : fi 4 


cosy 3 - COS ^ 0 ^ 

sin 2 0 o ^ 

2 

cos^ - cos 0 O ^ 
sin 2 0 o ' 


(28c) 


(28d) 


The coupling terms between the edge of first-order diffraction 


f 

the £ terms. For each pair of intersecting edges a different 
factor results due to these terms. Due to symmetry, one needs to 
consider the range 0° £ <f>' * 90° only for the plate rotation angle. 

This eliminates having to consider interactions between some edges. 
The remaining interactions that one must consider are. 

(1) lst-order diffraction from edge 1 to edge 4. 

(2) lst-order diffraction from edge 1 to edge 2. 

(3) lst-order diffraction from edge 2 to edge 4. 

(4) lst-order diffraction from edge 2 to edge 1. 

(5) lst-order diffraction from edge 3 to edge 2. 

(6) lst-order diffraction from edge 3 to edge 4. 

(7) lst-order diffraction from edge 4 to edge 2. 

(8) lst-order diffraction from edge 4 to edge 3. 

The |cr X t 2 | terms for each of these interactions reduce to: 

Edge 1 to Edge 4: X I = cos ^o^ (29a) 

Edge 1 to Edge 2: | X ^2 2 1 (29b) 

Edge 2 to Edge 4: |cr 2 X t 2 J = cos/3 0 ^ (29c] 
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Edge 

2 

to 

Edge 

1 : 

1*2 

X 

A 

t2 J 

= sin0 o ' 

(2 9d) 

Edge 

3 

to 

Edge 

2 : 

1*3 

X 

A 

t 2 2 l 

= cosp 0 ' 3 

(29e) 

Edge 

3 

to 

Edge 

4 : 

A 

1*3 

X 

A 

- sinfV 3 

(2 9f ) 

Edge 

4 

to 

Edge 

2: 

A 

1*4 

X 

A 

t 2 2 l 

= COS0 O ' 

(29g) 

Edge 

4 

to 

Edge 

3: 

A 

1*4 

X 

A 

t2 3 l 

= sin|3 0 ' 

(29h) 


The distance parameter, £, designates the distance from the point 
of first-order diffraction to the point of second-order diffraction 
measured along cr, the grazing diffracted ray. The £ parameters are 
constant functions of incidence angle only for opposite edge 
interactions but are functions of distance along the edge for adjacent 
edge interactions. Fig. 5 shows the geometry for determining the £ 
parameters for interactions between edges 1 and 2 and edges 1 and 4. 
The geometries for the other interactions are similar. The limits of 
integration in the vector potential integrals of (15) vary according 
to the illumination of the edge of second-order diffraction. The 
extent of illumination is bounded by the grazing diffracted rays from 
the edge of first-order diffraction. Thus, the limits of integration 
are a function of the incidence angle. The distance parameters, £ and 
L, along with the limits of integration are given below for each pair 


of interacting edges. Recall that: 
L - ki£ sin 2 /3 Q ' 


(30) 
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(a.b) 


Edge 1 

- X 



(31a) 


Edge 1 to Edge £ 


£ 


14 


a - x , 


cos0' 


1*24 — Poj 

Limits of integration: “ a - x “ a 

a - 2btan(n-|3 ' ) i x s a 

Edge 1_ to Edge 2_ 

£ = — — 

12 cos|3 0 ' i 

L^2 = ^^12 

Limits of integration: 


-b 5 y ^ b - 


2a 


tan 


Edge 2 to Edge 4_ 


a + x 4 
” sin^o 2 

2 

L»24 = kJ£ 14 sin f 3 


Limits of integration: 


-a s x ^ a 


for 0 O ' ^ 7T - tan (a/b) 
for 0 O ' £ n - tan 1 (a/b) 


for p Q ' ^ 7T - tan 1 (a/b) 


for |3 0 ' — tan (a/b) 


< X s -a + 2btan0 o ' for 0 O ' * it - tan (a/b) 


Edge 2 to Edge 1 


£ 


2a 


21 sinfi 0 ' 2 


Lji = k£ 2 l sin 2 0 o ' 


Limits of integration: 

-b - y - b - 


2a 

tan (0 O ' ) 


for 0 ' ^ tan' (a/b) 


(31b) 


(31c) 

(31d) 


( 31e) 
( 3 1 f ) 


(31g) 

(31h) 
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{ 3 1 i ) 


Edge 3 to Edge 2 


i£ 


b - y 2 

32 ” sin/3'~ 


L32 “ |3 

Limits of integration: _ b _ y — b 


b - 2 atan^ 0 ^ s y * b 


Edge 3 to Edge £ 


£34 


2 b 


“'34 


sin£ ' 


kie 34 sin 0 ' 


Limits of integration: 

2b 

tan/3' 


a - 


< x i a 


Edge 4 to^ Edge 2_ 


£*■> - 


b + y 2 


42 COS /3 o ; 

L42 = k^ 42 S ^ 11 3 o 4 

Limits of integration: -b i y ^ b 


for 0 O ' ^ tan 1 (b/a) 
for /3 q ' 3 i tan " 1 (b/a) 


-1 


for / 3 0 ' ^ tan" (b/a) 


-1 


for 0 O ' i ti - tan (b/a) 


-b i y i -b + 2 atan (n-f3 ') for 0 O ' 4 * tan " 1 (b/a) 


Edge _4 to Edge 2 


= 


2 b 


43 cos/3 0 ' 

L 43 = k ^ 43 sin 2 |3 Q ' 


Limits of integration: 

2 b 


tan (Tr -0 ' ) 


< x i a 


for 0 O ' - 71 ' tan 


-1 


(31k) 

(311) 


(31m) 

(31n) 


(31o) 

(31p) 

(b/a) 
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The second-order fields are determined by substituting the 
current components for each pair of interacting edges into the vector 
potential equations of (ID and using (10). For opposite edges the 
integrals reduce to closed-form expressions. The adjacent-edge 
integrals must be evaluated numerically. Just as for the first-order 
diffractions, integration is with respect to either the x o y 
coordinate so that the integrals involved simplify to the following: 

Edge 1_ to Edge _4 
_a = X j 


14 


14 


w jkxsinG' cos <P r , 

e " jk ^i4 F(/ L 14 (l-jiiY ) e dx 


-) 


..(ft jkxsinG' cos</>' , 

e f (V L 14 ( 1+cosi//^) ) e 


Edge 1 to Edge 2 


L 12 


_ ^ j2kysin0' sin <p' ^ 


Vl 


exp 


L, l±jl sin0'sin<£’l 

{ ^ t an (7t-/3 0 ' ) J 


= - 2 bsinc (2kbsin0' sin</>' ) - j2k sin0' sin#' 


(32a) 

(32b) 


(32c) 


Edge 2 to Edge 4. 


24 


24 


e "jk^24 f (V L 24 ( 1-M2> > e 


jkxsinG' cos<p r 


-1 

X 1 

x 2 . ■ jkxs inG' cos^' , 

= J e '^ kie 24 F (i/ L 24 (l+cosi// 2 ) ) e dx 


(32d) 


(32e) 
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Edge 2 to Edge 1 


21 


^ 12 kysin 0 'sin <b' 

J e dy 


Vl 


exp -j 


4ak 


= 2bsinc {2kbsin0' sin<£' ) + 


tan/3 o : 


si 


n 0 ' sin </>' j 


j2k sin0' sin<£' 


( 3 2 f ) 


Edge 3 to Edge 2 


32 


- '/ 


- jk jE 
e 

*1 


„ . r~ rr — . jkysin 0 ' sin 0 ' , 

32 F (v L 32 ( ) ) e J J dy 


y 2 


1° = -/ e~ jk ^ 32 F (V L 32 (I+COS 1 // 3 ) ) 

32 *i 


jkysin0'sin0' , 

e dy 


Edge 3 to Edge 4 


34 


X 9 

7 j 2 kxsin 0 ' cos©' 

J e dx 


Edge 4 to Edge 2 


f 4bk ^ ^ 1 

exp -j — s inS'cos^' I 

l tan p 0 ; J 


j2k sin0' cos<p' 


< 2 TiV ) e lkl ” 1 " 9 '- in ' > ' dy 


I* - -f 

42 / 

Vl 

I /j2 = -f e 3 42 F (V L 42 (l+cosi// 4 ) ) e 

*1 


jkysinO' sin<f>' 


dy 


Edge 4 to Edge 3 


L 3 4 


^ j 2 kxsin 0 ' sind' 

J e dx 


exp j 


4 kb 


tan (JT -0 ' ) 


- — sinG'cos^' 


j2k sin0' cos 0' 


( 32g ) 
( 32h) 


( 3 2 i ) 


(32 j) 
(32k) 


(321) 
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B. Corner Diffraction 


Because the GTD and UTD diffraction coefficients are derived from 
the exact solution to scattering by an infinite wedge, the 
coefficients fail to account for the joining of two edges at a corner. 
For certain aspect angles, the scattering from the corners is 
significant. Pathak and Burnside developed a heuristic corner 
diffraction coefficient [7] -[8] based upon an appropriate, although 
non-rigorous, modification of an asymptotic evaluation of the 
radiation integral due to the equivalent edge current that would exist 
along the scattering edge if the corner were not present. One major 
flaw in the coefficient is that it is non-unique for certain 
backscatter angles near normal incidence [9]/ thus a 
rigorously-developed corner diffraction coefficient is desirable. 
However, Pathak and Burnside's coefficient is successful for many 
plate geometries and may be used with caution. 

The geometry for a corner in a planar surface is shown in Fig. 6. 
The total diffracted field from one corner is the sum of contributions 
from each of the edges comprising the corner. The general form of the 
corner diffracted field is: 


E C (s)=-E 1 (Q c ) -D c - /- 

— — ^ L v s (s +s ) V 


S (S+S c ) 


- j ks 


(33) 


where 

E 1 (Q c ) - the incident field at the corner. 

D c = the dyadic corner diffraction coefficient. 
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s = the distance from the corner to the receiver. 

s' = the distance from the source to the point of edge 
diffraction . 

s" = the distance from the point of edge diffraction to the 
receiver . 

s c = the distance from the source to the corner. 

“ k, s 

e * the phase factor. 

The dyadic corner diffraction coefficient, like the ordinary 
diffraction coefficient, is in terms of parallel and perpendicular 


components : 


A i A dc c A i A dc c 
D c = e h e h D h + e s e s D s 


(34) 


and Dg are the hard and soft corner diffraction coefficients. 


respectively, given by the following: 
- j7T/4 


D 


s , h 


V 2 TTk 


h (Qp.) 


v /sin(3 0 sin|3 0< 


s,h '“ E ' co S0 OC - COS^ c 


F [kL a (rr+|3 - 0 C > ] 


(35) 


where 


C s,h 


- j 71/ 4 

-e 

2V 2uk sin/3 0 


x 


r 

F [kLa O') ] 

r La O ) /X 1 

cos O / 2 ) 

[kL c a (ir+/3 oc - P C )J 


_ F [kLa 0 + ) ] 

r La (/3 + ) /A 1 


+ 

cos 0 + /2) 

[kL c a(n+p oc - 0 C )J 

J 


( 36 ) 
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F is the Fresnel transition 


function given by: 


, jx f -jT 

F(x) = 2yf x e J e 

VZ 


dr 


( 37 ) 


j3 are as: 

+ (38) 

(3 = ip ± «//' 

0 o is the Keller cone angle. The other angles and functions are: 

p c = the angle between the incident ray at the corner and 
the edae of interest. 

R = the angle between the diffracted ray at the corner 

° and the extension of the edge of interest as s ow 

Fig. 3-5. 


2 

a (\p) - 2cos ( i/// 2 ) 


L = 


S' S’ 


(s' +s") 


sin (|3 0 ) 


s c s 

L c “ S c +S 


(39a) 

(39b) 

(39c) 


These fields simplify considerably for far field, plane wave 

scattering from the rectangular plate. 

The RCS of the flat plate in all planes can be determined using 

only the corner diffraction coefficient. Near and at normal 

incidence, this formulation fails due to the nonreciprocal nature of 
the corner diffraction coefficient. To alleviate this problem, near 
and at normal incidence, the GTD equivalent currents solution is used; 
and the corner diffraction results are used away from the problem 
area. The total scattered field from a corner consists of terms for 
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each of the adjoining edges; therefore, eight terms are needed for the 
rectangular plate . 

The general expression for the corner diffracted field due to one 
of the two adjoining edges is given in (33) . For far-field 
scattering, s e s c = s' = s" ^ » ; therefore: 

- j ks 


E C (s) = -E 1 (Q c ) • D c 


for far-field scattering. 

For far-field backscatter ing, L c = oo so that: 
F[kL c a(7t+p oc -p c ) ] 8f 1 
Also, /3 oc = n - (3 C and ip' = \p so that: 


(41) 


(42) 


V sin/3 sin/3 


oc -tan/3 


c °s (3 oc - cos /3_ 


(43) 


The diffraction coefficient simplifies to: 

- jTT/4 


D 


-e 


s , h 


C.w(Q tt ) 


tan/3 


, h "•£' — 2 


(44) 


V 2nk 

For far-field backscattering, L = L c = oo and C sh (Q E ) simplifies to: 

- j n/ 4 


C s ,h(Q E ) = 


2v2Trk sin/3 c 


{—+ ) 

^ 2TTCOS (7T“P 0 ) ' 

‘ u i — I 

.osip | [ 2ttcos 2 (n-/3 Q ) ' J 


(45) 


COS x}) 

The total field scattered due to one edge adjoining at a corner 
simplifies to: 
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E c s , h < s > - E <Q c> 87ikc8¥p: 


i F r i .i. 

\ 2 _ I COS 

[ ' 2ttcos /3 0 } 


cos 0 
2 

2ncos f3 ( 


The top view of the plate geometry for the corner diffraction 
analysis is shown in Fig. 7. The incident fields are the same as in 
(9). The field must be determined at the corner of interest. The 
corners are designated as the following: 


Corner A: 


Corner B: 


Corner C: 


Corner D: 


x “ a 


x = -a 


x = -a 


x = a 


y *= b 
y = b 
y = -b 


The angles 0 are the same as those given in (5) . The * angle used for 
the scattering component for a designated edge and corner is the «/» 
associated with the edge. For example, the angle is used for 

scattering from corner A due to the presence of edge 1. 

A 

The |3 0 angle is the angle between the -s' and the edge or: 


cos/3 0 “ -s ' • c 


where 


c = the unit vector tangent to the edge of interest pointing 
outward from the corner of interest. 

The vectors and angular functions for each corner/edge combination 
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Corner A 


Edge 1: c A1 - 

A 

” a y 

c °S^ 0fil = - 

sin0' sin#' 

(48a) 

A 

Edge 3 : c A3 = 

A 

~ a x 

COS ^A3 = ' 

■sin0' cos#' 

(48b) 

Corner B 





Edge 2 : c B2 ” 

A 

-a y 

COSl3 °B2 = ' 

-sin0' sin#' 

(48c) 

A 

Edge 3 : c B3 = 

/A 

a x 

c °sPo B 3 = 

sin0' cos#' 

(48d) 

Corner C 





A 

Edge 2 : c c2 =* 

A 

a y 

COS ' 3 o c2 = 

sin0' sin#' 

(48e) 

A 

Edge 4 : c c4 - 

A 

a x 

cos Po c ^ 

sinO' cos#' 

(48f ) 

Corner D 





A 

Edge 1 : c D1 - 

/\ 

a y 

c °s/ 3 0d1 “ 

sinO' sin#' 

(48g) 

A 

Edge 4 : c D4 = 

A 

a x 

cos Po D4 = 

-sin0' cos#' 

( 4 8h) 

The final parameter that 

must 

be designated is s, the 

distance from 

the corner to the observation 

point. The 

geometry indicating these 

distances is in Fig. 8. 

For 

far-field scattering, the 

following 


approximations are used: 


Amplitude : 

• S A " S B “ S C 

- Sp, - s 


(49a) 

Phase : 

sin0' / 

S A - s 2 / 

2(a 2 +b 2 ) 

(cos#' +sin#' ) 

(49b) 


sin0' / 
S B ~ 3 - 2 / 

2 (a 2 +b 2 ) 

(sin#' -cos#' ) 

(49c) 


sin0' / 

s c “ s + 2 v 

2 (a 2 +b 2 ) 

(cos#' +sin#' ) 

( 4 9d) 


sin0' I 

S D - s 2 / 

2 (a 2 +b 2 ) 

(cos#' -sin#' ) 

( 4 9f ) 
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The total field scattered by the plate consists of 
of (46) for each edge joined at a corner. The=e 
added to arrive at the total field. 


a term 
eight 


of the form 
terms are 
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III. Results 


Computations were made for a square plate with each side equal to 
5.73 A. Comparisons are made among the following: mm and 

experimental results, the PO equivalent currents model, the GTD 
equivalent currents model, the PO/PTD equivalent currents model, the 

PO/PTD/ 2nd- order equivalent currents model, and the corner diffraction 
coefficients model. 

Fig. 9 shows soft polarization results for a 30° rotated plate. 
Even the PO equivalent currents, which account only for surface 
scattering, give good results near normal incidence. As the angle of 
incidence moves away from normal, there is a need for components to 
account for edge diffraction. The GTD and PO/PTD models greatly 
improve the results in the grazing regions, although there remains 
some disagreement which points to the necessity for higher-order and 
corner diffraction components. Fig. 10 shows results from the PO/PTD 
model with the second-order PTD coefficients added and also shows 
results from the corner diffraction coefficient model. Since the 
diffraction coefficients are inaccurate near and at normal 
incidence, the GTD equivalent currents solution is used in the region 
±5° on either side of the normal direction and the corner diffraction 
solution is used elsewhere. The second-order currents do not improve 
the results for this case and even result in worse agreement than that 
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obtained using the first-order models. The corner diffraction model 
improves the results near grating and agrees »ith the first-order 
models away from grating, indicating that the corner diffraction 
mechanism is the more crucial scattering mechanism for this plate 
rotation and polarization. 

Fig. 11 shows the hard polarization results using the first-order 
models for the same plate configuration rotated 30°. Near and at 
grazing there is a major discrepancy between the first-order, 
high-frequency models and the experimental and MM results. The PO/PTD 
models with the second-order components added and the corner 
diffraction coefficients model yield much better results although 
discrepancies still exist. These results are displayed in Fig. 12. 

O 

The soft polarization results for the same size plate rotated 45 
are in Fig. 13. One would expect that corner-diffraction would play a 
major role at this angle of rotation so that the large discrepancies 
of the first-order models near grazing incidence are not surp g 

The addition of second-order equivalent currents yields excellent 
results in this case. Surprisingly, the corner diffraction model does 
not improve the results. These results, shown in Fig. 14, indicate 
that for this angle of rotation and polarization second-order 

components are the major contributing factors. 

Hard plarization results for the same rotation angle are shown in 
Figs. 15 and 16. Just as for the 30° rotated plate for hard 
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polarization, the first-order, high-frequency models do not produce 
good results near the grazing regions. The corner diffraction model 
greatly improves the results; however, the addition of second order 
components does not result in an improved model. This indicates that 
the corner diffraction mechanism is dominant for this configuration. 

IV. Conclusions 


The nonprincipal-plane scattering from a rectangular flat plate 
was considered. Comparisons among five high-frequency models, MM and 
experimental results were made. Near normal incidence all the models 
agreed; however, near grazing incidence a need for higher-order 
and corner diffraction mechanisms was noted. In many instances the 
second-order and corner-scattered fields formulated in this report 
improved the results. 
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PATTERN CONTROL OF HORN ANTENNAS 
I. Tasks Accomplished 

During this period, the computations of the impedance elements 
have been completed. These include interactions between the two 
electric current modes, the electric current mode and the magnetic 
current mode, and the two magnetic current modes. Especially, an 
accurate and efficient formulation of computing interactions between 
electric current mode and magnetic current mode has been accomplished. 
This, together with other subroutines we have developed, enables us to 
fill-in all the elememt in the matrix. 

After the fill-in of the impedance elements in the matrix, the 
forward problem is accomplished. That is, given the specification of 
the horn and the excitating waveguide mode, the radiation pattern of 
the antenna based on the integral equation can be obtained. 

An example case was run for a standard X-band gain-horn (DBG-520) 
with the configuration in Figure 1. The H- and E-plane patterns of 
this horn antenna with perfectly conducting walls are shown 
respectively in Figures 2 and 3. Comparison with the gain pattern 
available from the manufacturer for up to the first side lobe shows 
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Figure 2. Normalized amplitude pattern for 20dB gain horn 
antenna ( DBG — 520 ) with perfectly conducting electric 
walls ( f = 10 GHz ). 


Plane Pattern 
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Figure 3. Normalized amplitude pattern for 20dB gain horn 
antenna ( DBG — 520 ) with perfectly conducting electric 
walls ( f - 10 GHz ). 


QF foor Quality 


good agreements, although the cross polarization has not yet been 
accounted for . 

To investigate the effect of the lossy coating on the radiation 
pattern, two sets of the lossy materials were used to cover the top 
and bottom walls of the horn. This is intended to improve the E-plane 
pattern. The plots shown in Figures 4 and 5 are obtained by uniformly 
covering the top and bottom walls with a layer of AlSb with a 
thickness of 0.001A o ( 3 x 10 ^meters) . The material has a relative 
dielectric constant of 11 and a resistivity of 0.005 C2*m. The 
resulting E-plane pattern shows about 2-dB improvement in the 
sidelobes. Figure 6 and Figure 7 are obtained based on a material 
which has a relative dielectric constant of 3 and a sheet resistance 
of 1500 Q per square. The thickness of this material is 2 mils (5.08 
x 10 ^ meters ) . The resulting E-plane pattern shows about 3-dB 
improvement in the first sidelobe and 4-dB improvement in the second 
sidelobe . 

Further improvement can be expected by increasing the thickness 
of the coating. However, the validity of our impedance boundary 
condition becomes questionable. A better impedance condition has been 
developed; however it has not yet been implemented in the computer 
program. 

In the moment method solution of this project, the computations 
of the matrix elements are the most tedious part of the work. We have 
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Figure 4. Normalized amplitude pattern for 20dB gain— horn 
antenna ( DBG — 520 ) with upper and lower walls covered with 
AlSb ( f = 10 GHz). 
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Figure 5. Normalized amplitude pattern for 20dB gain— horn 

an ^ enna ( DBG 520 ) with upper and lower walls covered with 
AlSb ( f = 10 GHz). 
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Figure 6. Normalized amplitude pattern for 20dB gain— horn 
antenna ( DBG— 520 ) with upper and lower walls covered with 
lossy material ( f = 10 GHz). 
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basically fullfilled this task. Although we still did not have time 
to connect all parts of the work together to realize the synthesis 
problem, we can say we are progressing well toward that goal. 


II. Future Work 


Future work will be concentrated on the following items: 

To include the cross polarization components of the 
equivalent magnetic currents on the two appertures . 

To compare this integral equation method with another 
rigorous method by H. Patzelt and F. Arndt [1], 

J To investigate the realization of the sheet impedance needed 
to control the radiation pattern, and to extensively verify 
the validity of the impedance boundary condition. 

4. To reassemble the matrix equation to solve the synthesis 
problem as was presented in the previous report. 
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